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Abstract 


A dynamic one-dimensional rigorous process model of a single-cell direct methanol fuel cell (DMFC) is presented. Multi-component mass 
transport in the diffusion layers and the polymer electrolyte membrane (PEM) is described using the generalised Maxwell—Stefan equation for 
porous structures. In the PEM, local swelling behaviour and non-idealities are accounted for by a Flory—Huggins activity model. This model is 
used as basis of a model family with different anode and cathode reaction mechanisms (single-step and multi-step with and without adsorption to 
catalyst surface sites). 

The model variants were used to simulate the dynamic (transient) response of the DMFC to stepwise changes in the methanol feed concentration 
from typical operating levels down to zero, while maintaining the cell current. For validation, similar experiments were carried out. In the 
experiments, the cell voltage broke down only after an unexpectedly long period of time, and for a variety of operating conditions even a cell 
voltage overshoot could be observed. Such overshoot behaviour is also predicted by those model variants, which feature anode reaction mechanisms 
with reaction intermediates (e.g. CO) adsorbed to the anode catalyst, while models without such detailed anode reaction mechanisms fail in this 
respect. The model-based analysis reveals that the observed overshoots result from the different time constants of the responses of the anode and 


cathode overpotentials to the feed change. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The direct methanol fuel cell (DMFC) emerges as an interest- 
ing technology especially for portable power supply units. For 
such portable systems, dynamic rather than stationary operation 
is typical. But not only possibly rapid changes in the electrical 
power demand have to be considered, but also the fact that such 
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systems have to be quite simple, which, e.g. makes it difficult to 
keep the methanol feed concentration constant. Cheap and reali- 
able sensors for methanol concentration are not available, and 
therefore a DMFC system will have to cope with changes in the 
methanol feed concentration, besides other operating conditions. 
Moreover, the examination of transient responses of a system to 
changes in the operating conditions can yield deeper insight into 
the governing internal physico-chemical phenomena. 

In this work, the response of a DMFC to step changes in the 
methanol feed concentration is examined. Systematic experi- 
ments revealed a very interesting dynamic behaviour especially 
when the anode feed is changed from methanol solutions to 
pure water while maintaining the cell current. For low current 
densities, even intermediary increases in the cell voltage can be 
observed (overshoot behaviour) [1,2]. 

In order to understand the governing physico-chemical phe- 
nomena, a rigorous dynamic model of a DMFC was formulated. 
To find out the role of the anode and cathode reaction mech- 
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Nomenclature 

a activity 

айо water vapour activity 

AS cell cross-sectional area (m?) 

B transport matrix 

с molar concentration in fluid phase (mol m~?) 

č molar pseudo-concentration w.r.t. total volume (in 
porous structures only) (mol m7?) 

d thickness, diameter (m) 

D Maxwell-Stefan binary diffusion coefficient 
(т? 5—1) 

е enthalpy flux density w.r.t. cross-sectional area Aê 
Jm? s7!) 

E^ activation energy (J mol~!) 

F Faraday’s constant, F=96,485 As mol`! 
(А s mol`!) 

h specific enthalpy (J mol~!) 

i current density w.r.t. cross-sectional area AS 
(А м7?) 

j individual molar flux density w.r.t. cross-sectional 
area AS (mol m~? s7!) 

k index for control volumes (discretised model) 

Li friction terms (s m7?) 

m mass flux density w.r.t. cross-sectional area AS 
(kg m~? 5—1) 

M molar mass (kg што!) 

п overall molar flux density w.r.t. cross-sectional 


area AS (mol m2 s7!) 

N number of moles (mol) 

NM,cu number of chain units between two polymer 
cross-links 


N mole density (loading) w.r.t. cross-sectional area 
AS (mol m~?) 

Np mole density of active Pt sites (anode/cathode cat- 
alyst) (mol m~?) 

p pressure (Pa) 

Psat saturation pressure (Pa) 

q heat flux density (due to thermal conduction) w.r.t. 
cross-sectional area AS (J m~? s7!) 

r reaction rate w.r.t. total pore volume 


(mol m7? s7!) 

ideal gas constant, R=8.314 J mol`! K7!) 
time (s) 

temperature (K) 

voltage (V) 

velocity (m s7!) 

volume (т?) 

molar volume (m? mol” !) 

mole fraction in liquid phase 

mole fraction in gas phase 

cell coordinate perpendicular to cell plane 
number of transferred electrons/single charges 


<<. GN7 


NNK X 


* 


Greek symbols 

Qa,% Charge transfer coefficients (anodic, cathodic) 

E€ volume fraction (pore volume fraction = porosity) 

n overpotential (V) 

МЫ dynamic viscosity (Pas) 

N thermal conductivity coefficient (W m`! K7!) 

A relative water content in membrane 

u chemical potential (J mol~!) 

v stoichiometric coefficient 

ф electrical potential (V) 

X non-ideality coefficient in Flory—Huggins activity 
model 

Superscripts 

A anode compartment (supply channel structure) 

AC anode catalyst layer 

ACP polymer phase within (AC) 

AD anode diffusion layer 

AF anode feed 

С cathode compartment (supply channel structure) 

CC cathode catalyst layer 

CCP polymer phase within (CC) 

CD cathode diffusion layer 

CF cathode feed 

eff effective 

M membrane (PEM) 

ө at standard conditions: T? = 298 К, p? = 105 Pa 

Subscripts 

a anode 

с cathode 

carbon carbon material 

cell cell 

cu polymer chain unit 


CH30H methanol 


СО» 
eff 
Ht 


carbon dioxide 
effective 
proton 

water 
counting index 
counting index 
Joule heating 
nitrogen 
oxygen 

in pore(s) 
saturated 
sound 


anisms in the overall DMFC behaviour, assumed reaction 
mechanisms with different complexity were implemented in a 
base model. The base model features a very detailed descrip- 
tion of mass transport phenomena, and has proven to be able 
to accurately predict steady state methanol and water crossover 
fluxes and cell voltages for a broad range of operating conditions 


[1,3,4]. 
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Comparison of the simulation results obtained from the 
model variants reveals hints on the role of reaction interme- 
diates on the platinum catalysts and methanol mass transport. It 
turns out, that intermediates adsorbed to the anode catalyst seem 
to play an important role, therefore an adequate dynamic model 
has to account for such phenomena. While similar mechanisms 
can be expected on the cathode catalyst as well (and are imple- 
mented in the most complex model variant), they seem to be of 
lower importance. 


2. Dynamic experiments 


Experiments were carried out using a single cell DMFC 
fed with air and liquid methanol water solutions. A detailed 
description of the DMFC design can be found in [1]. The 
identical anode and cathode monopolar plates are made from 
graphite material (thickness 7mm, material code FU4369 by 
Schunk Kohlenstofftechnik, Germany). The flowbed structures 
consist of parallel channels of 2mm width and 2mm depth, 
with 1 mm wide ribs between them. The flowbed itself has the 
outer dimensions 65 mm x 40 mm (active area of AS = 26 cm”). 
As diffusion layers PTFE-coated TORAY carbon paper (TGP- 
H-060) is used, with a PTFE loading between 20 and 25 mass% 
with respect to the uncoated material. The membrane electrode 
assemblies (MEA) are prepared from NAFION™ N-105 mem- 
brane foil, onto which the catalyst layers are applied using an 
airbrush technique. The anode catalyst layer features a cata- 
lyst loading of 5 тест? (unsupported) platinum ruthenium 
black (Alfa Aesar Johnson Matthey HiSPEC™ 6000) and a 
NAFION™ content of 15 mass% relative to the metal loading 
(i.e. 0.75 mg cm~?). The cathode catalyst layer has the same 
metal loading, but as catalyst (unsupported) platinum black is 
used (Alfa Aesar Johnson Matthey HiSPEC™ 1000) and the 
NAFION™ content is 10 mass% relative to the metal loading 
(i.e. 0.5 mg cm”). 

In the experiments, the methanol feed concentration of the 
cell was changed stepwise from different starting values down 
to zero (pure water feed), while all other operating parameters 
(temperature, pressures, flow rates, etc.) of the cell as well as 
the cell current were kept constant. Some typical results of these 
experiments are presented in Fig. 1. 
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Fig. 1. Methanol feed concentration step-down experiments (galvanostatic oper- 
ation) at different cell current densities and cell temperatures. 


The full set of operating conditions is: 


Anode 
Methanol feed concentration: 0н = 1/1.5 mol dm~? 
Pressure: p^ = 1.7 bara = 1.7 х 10° Pa 
Feed flow rate: F^ = 0.5 dm? min! 
Feed temperature: T^ = 60/75/90 °C 
Cathode 
Dry air (dew point approximately —2 °C) 
Pressure: p? = 1.7 bara = 1.7 х 10° Pa 
Feed flow rate: ЕСЕ = 0.5 scbmh7! (at 1.013 х 10° Pa and 20°C) 
Feed temperature: 7С = 30°C 


As the feed switch is not located directly at the fuel cell inlet, 
a certain dead time occurs between the moment the valves are 
operated and the moment the changed feed solution reaches the 
anode inlet of the DMFC. This dead time depends on the anode 
feed flow rate. The distance between the feed switch valves and 
the anode inlet is approximately 6m, with the pipes having an 
inner diameter of 6 mm. Assuming plug flow, this leads to a total 
dead volume of: 


л 
2 
Vdead = ApipeLpipe = д ripe pipe 


л 2 3 
= yom) 600 cm = 170 ст”. 


At the given flow rate of 500 cm? min™! the dead time results 
as: 


Vaead 170ст? 
FAF ~~ 500.cm3 min! 


ldead = = 20s. 


In Fig. 1, the dead time is marked by a vertical dashed line at 20 s 
after the feed switch. In all figures comparing experimental and 
simulation results, the dead time of the experiments is accounted 
for by substraction of 20s from real time since switch, as in the 
model no such dead time is accounted for. 

The dynamic response of the DMFC to step-down of the 
methanol feed concentration (Fig. 1) has two noteworthy fea- 
tures: 

First, the cell voltage does not break down immediately after 
the switch to pure water feed, but only after some 10 s, depending 
on operating conditions. This is unexpected, as the feed flow rate 
is very high in the experiments, so that the mean residence time 
of the anode liquid in the anode flowbed compartment is below 
1s. Therefore one can assume, that the anode compartment as 
the supposedly main reservoir for fresh methanol is emptied 
completely within less than a second. Obviously there must be 
other reservoirs for reactants within the inner structures of the 
DMFC [1]. 

Second, for low current densities the inevitable cell volt- 
age breakdown is not only taking place later than expected, 
but the cell voltage even shows an interim increase (overshoot 
behaviour) for low current densities (<100mAcm~?). Such 
experimental behaviour was already reported elsewhere [2]. 

The results of the systematic experiments presented in 
Figs. 1-3 can be summarized as follows: 
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Fig. 2. Methanol feed concentration step-down experiments (galvanostatic oper- 
ation) at different cell temperatures. 


e In all cases (with or without overshooting), the cell breaks 
down much slower than the anode compartment is rinsed from 
methanol. 

e The sustaining time fsustain (i.e. the time between end of dead 
time fdeaq and the time the cell voltage has broken down to 
zero) is only a function of the cell current density, but not of 
the methanol feed concentration and cell temperature. 

e Overshooting of the cell voltage only occurs for low current 
densities (Fig. 1). 

e The higher the cell temperature, the higher the level of over- 
shooting (Fig. 2). 

e The higher the initial methanol feed concentration, the higher 
the level of overshooting (Fig. 3). 
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Fig. 3. Methanol feed concentration step-down experiments (galvanostatic oper- 
ation) at different methanol feed concentrations and cell temperatures. 


Explanations for this behaviour are based on the combina- 
tion of several assumed physico-chemical phenomena occuring 
simultaneously in the DMFC: 


e Methanol inside the pores of the anode diffusion layer and 
anode catalyst layer is not immediately rinsed out of the 
DMEC, therefore the methanol concentration in the anode cat- 
alyst layer is not decreasing as fast as that in the anode flowbed 
structure. As already shown elsewhere [1], this amount of 
methanol is too small to explain the observed sustaining times. 

e Methanol inside the PEM could reverse its flow from cathode 
towards anode when the anode is completely rinsed, and thus 
supply the anode reaction with additional reactant. 

e Reaction intermediates (e.g. CO) adsorbed to the anode cat- 
alyst are another possible reactant reservoir. If the supply 
of fresh methanol ceases, such intermediates could be fully 
oxidised, thus allowing further electric current flow. 

e Decreased methanol crossover, due to decreasing methanol 
concentration in the anode catalyst layer, would lead to a 
decrease in the absolute cathode overpotential. 


3. Model formulation 


To analyse the experimental behaviour of the DMFC, a one- 
dimensional rigorous dynamic process model of a DMFC was 
developed [1,3,4]. It reflects the seven-layer structure of the 
DMFC: anode compartment (A), anode diffusion layer (AD), 
anode catalyst layer (AC), polymer electrolyte membrane (M), 
cathode catalyst layer (CC), cathode diffusion layer (CD) and 
cathode compartment (С). For all these layers, dynamic mass and 
energy balances, as well as quasi-steady state charge balances for 
(AC) and (CC) are formulated, and appropriate transport kinet- 
ics and internal boundary conditions are used to couple them. 
The diffusion layers and the PEM are spatially distributed ele- 
ments, while the other four elements (catalyst layers and supply 
compartments) are described as ideally mixed systems. Fig. 4 
presents a scheme of the model. The spatially distributed ele- 
ments are described along the coordinates zAD zM апа z®, 
respectively, where the coordinates are perpendicular to the 
cell plane and increase from anode to cathode side. The model 
has proven to be capable of predicting measured steady state 
methanol and water crossover fluxes through the PEM under var- 
ious operating conditions, with only one set of parameters [1,3]. 
For low to moderate current densities also the steady state current 
voltage characteristics are predicted with acceptable accuracy. 
A brief overview of all governing equations (balances, transport 
kinetics) is given in Table 1 (anode and cathode compartments, 
diffusion layers and catalyst layers) and Table 2 (PEM). The 
equations connected with the anode and cathode reaction kinet- 
ics (Table 3) will be discussed in more detail in the following 
sections. The symbols used in Tables 1—3 are explained in the 
list of symbols. 

In this contribution, the model is used to simulate the dynamic 
response of the DMFC to changes in the methanol feed concen- 
tration. The analysis focuses on the role of the anode and cathode 
reaction kinetics, for which purpose several different reaction 
kinetic approaches for anode and cathode (Table 3) were imple- 
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Fig. 4. Working pronciple of DMFC, model structure and geometry. 


mented in the base model (Tables 1 and 2). The simulation results 
of these different model variants are then compared to experi- 
mental data in order to find out, which influence the electrode 
kinetics have on the overall dynamic behaviour of the DMFC. 
The final aim is to identify appropriate anode and cathode kinetic 
models, which can predict and physico-chemically explain the 
observed dynamic behaviour. 


3.1. Model I: base model 


In the base model (in the following referred to as model I), 
lumped reaction schemes are assumed [1,3]. Table 3 presents 
those reaction schemes (anode: T3.1; cathode: T3.12 and T3.13), 
the mass and charge balances (anode: T3.2 and T3.3; cath- 
ode: T3.14 and T3.15), and the corresponding rate expressions 
(anode: T3.4; cathode: T3.16 and T3.17). 

The anodic methanol oxidation (Eq. T3.1) as well as 
the cathodic oxygen reduction (Eq. T3.12) are assumed to 
be reversible single-step electrochemical reactions, conse- 
quently their kinetics are decribed by classical Butler—Volmer 
approaches (Eqs. T3.4 and T3.16). At the cathode a second, 
undesired reaction takes place: methanol permeating through the 
PEM is oxidised to carbon dioxide (Eq. T3.13). It is assumed 
that this reaction is very fast, so that all methanol reaching the 
cathode is immediately consumed, independent of the electrode 
potential. Therefore, the rate of this reaction is proportional to 
the methanol flux reaching the cathode (Eq. T3.17). 

The two reactions at the cathode are coupled via the com- 
ponent mass balances (Eq. T3.14) and also via the quasi-steady 
state charge balance (Eq. T3.15). The oxygen reduction con- 
sumes protons and electrons, while the methanol oxidation 
produces them, therefore the electric potential of the cathode 
is a mixed potential of both reactions’ contributions. 


3.2. Model IT: extended anode kinetics 


It is well known from DMFC literature that the anodic 
methanol oxidation on platinum ruthenium catalyst is neither 


a single-step nor a reversible reaction. It consists of various par- 
allel and consecutive reaction and adsorption/desorption steps, 
with quite stable reaction intermediates bound to the different 
catalyst surfaces (usually platinum and ruthenium) [5]. Still the 
reaction mechanism is not fully understood and various possi- 
ble mechanisms and corresponding rate expressions have been 
proposed, e.g.[5—8]. 

To account for the multi-step nature of the anodic methanol 
oxidation while limiting the complexity of the reaction scheme, 
a two-step mechanism is used here, which is a slightly sim- 
plified version of one that has been published recently [9]. It 
is presented in Table 3 (lower left quarter). In a first, irre- 
versible step, methanol is adsorbed electrochemically to the 
platinum catalyst, forming strongly adsorbed carbon monox- 
ide, four protons and four electrons (Eq. T3.5). The second step 
(T3.6) is a combination of three consecutive reactions: (a) water 
adsorbes electrochemically to the ruthenium catalyst in a fast 
and reversible reaction. (b) The formed Ru—OH then reacts irre- 
versibly with the CO adsorbed on the adjacent platinum sites 
to carbon dioxide, two protons and two electrons. (c) Finally, 
the formed carbon dioxide desorbs from the electrode surface. 
In these three reaction steps, the irreversible reaction between 
adsorbed water and adsorbed CO (b) is assumed to be the rate 
determining step, while the water adsorption on ruthenium (a) 
is assumed to be always in quasi-equilibrium, and the carbon 
dioxide desorption (c) is assumed to be a fast and irreversible 
reaction. 

Resulting from this reaction scheme, an additional compo- 
nent mass balance for the surface coverage of the platinum 
catalyst with CO (OS, Eq. T3.8) is necessary, and the charge 
balance (T3.9) has to account for both reactions’ charge produc- 
tions. 

As both combined reaction steps are irreversible electrochem- 
ical reactions, their rate expressions (Eqs. T3.10 and T3.11) 
are formulated as Tafel kinetics. Both reactions also include 
adsorption and desorption of species to the platinum cata- 
lyst surface, so the rate expressions also contain a term to 


Table 1 


Model equations and parameters of base model (flowbed compartments, diffusion layers and catalyst layers) 


Components Anode compartment (A) 


Index j= НО, CH30H, СО» (pure liquid phase) 


Cathode compartment (C) 
Index j = №, O2, H20, СО» (pure gas phase) 


dc 1 dc 1 
Component mass balance a VA (BAF (АВ ch) | Asn}? | Ар} (T1.1) a үе (FES ЕС } А9 | .cp_ycp) (ТІ.13) 
ah ЕС cr, КТС s cD 
with F~“ = Е + СА п} [Ке гез] (Т1.14) 
Ј 
Total mass balance No balance, pressure р^ is operating parameter! No balance, pressure рС is operating parameter! 
Energy balance No balance, temperature T^ is operating parameter! No balance, temperatures 7< and ТС carbon ате Operating parameters! 
Components Anode diffusion layer (AD) Cathode diffusion layer (CD) 
Index j=H20, CH30H, СО» (pure liquid phase) Index j = №, O2, H20, СО» (pure gas phase) 
„Ар AD .АР__ AD „Ар ‚АР ac 1 дпСР 
Component mass balance Ар д Xj =>. Mi dj х) Si ws (ТІ.2) 5 Ср z (TL15) 
tot z TT эг? ‚еў Din eff Epores zZ 
a“? _ 2 = 
Total mass balance ЕН —(v4 Veoun D о (Т1.3) No balance, ideal gas 
2 рб = КТ X с (TL16) 
apap 1 jeAD — ggAD әт©Р 1 Je деСР 
Energy balance - ы pa (TL4) - £ Ч (ТІ.17) 
дг ~ AD д2 az т az д2 
(pcp) (pcp) 
AD AD AD 7-1 AD CD 1 Ср ү-1 CD 
)=-|У c IBPT V (x) (=e lB T V (p°) 
т КТ 
Individual component flux densities Б Bi? 1 AD Re? xf? (ТІ.5) Т ср 1 Ў “Dees B? у? (Т1.18) 
(Maxwell-Stefan mass transport model) Wit ogy ар > jir) ADF wi “Ж Sor = Во 
Bi poet fi pee pee DZE A pas. в 
у ) apa (а)? ap 
Convective flow velocity уАр = i R (TL6) v? = Laks р (Т1.19) 
P КШ dz р 19 д: 
Total component flux densities nv = = ЈА б Dis тоо with ae = с^, бле (ТІ.7) п = = с р D + oe ур with C= 2а (T1.20) 
Total mass flux density mad = =2 (n}? Mj ) (T1.8) Not necessary as no total mass balance formulated 
Total enthalpy flux density р = J (п^0 л T) (ТІ.9) = -2 п; (TO) (T1.21) 
AD Ар,ет TAP CD _ _4CD.eff are? 
Conductive heat flow density ge = А Ета (ТІ.10) q ле? Pa (T122) 
Z Z 
Components Anode catalyst layer (AC) Cathode catalyst layer (CC) 
Index j = НО, CH30H, СО» (pure liquid phase) Index j =N2, O2, H20, СО» (pure gas phase) 
Component mass balance See individual model formulations! See individual model formulations! 
dpac v2 
Total mass balance L 9 (тм, Meo mAb |2 AD _gAD ) (ТІ.11) No balance, ideal gas law: 
pee = RTS ese (T1.23) 
j 
dTAC (Ml Mo — e^Pl angan) + (lM = 2? | ap_ gab) атсс (e| coo — MI mgm) + (Pl coo = 2M Ima) 
Energy balance a cat 20 те cane а (ТІ.12) а; лш - = СС 20 cities (ТІ.24) 
(осу) а^ (осу) 4% 
Charge balance See individual model formulations! See individual model formulations! 


Rate expressions See individual model formulations! 


See individual model formulations! 
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Table 2 


Model equations and parameters of base model (polymer electrolyte membrane) 


Components 


Polymer electrolyte membrane (M) 
Index j= H,0, CH3OH, (Н?) (pure liquid phase) 


Component mass balance 


Energy balance 


Charge balance 
Total component flux densities 


(Maxwell-Stefan mass 
transport model) 


Phase equilibrium (AC-ACP) 


айу, м ; N Nj 
Тт = п рд with №; = АЗ (Т2.1) 

атм 1 дем ƏM agM 

(осу) Se í 

aiM 

0= – a (T2.3) 
M M ісер 
анъ = р F (T2.4) 


N Ін,о + Lengon — пн+ Li — nm oL2 
"CH30H Ls (T2.5) 


М Lmo = пн+ (La — 1016/13) — (L2 + L3)Lo/L3 


n (Т2.6) 
H20 
2 Ls — 1516/13 
о ён. Chon CHOH 
: _ o “но га; 3 з 
with Lino = M a > Loon = — = ЭЕ СУ 
H20 CH3 OH 
Хп,о Хеп,оп х 1 
= 2 ` Xir 
p eS А 
pI Э? 2 pt Э? 
H*, H,O H*,CH,OH , нҥ.но “ном, 
Хн 1 Хһ,о 
L=——_ + ——___ L, = - A 
3 pt p! 4 pt 
H*,CH,OH CH,OH,M , H',H,0 
Tee х Хоп,оп 1 L Xn,0 
б pe pf 6 per 
Brno Эң о,сн,он Pio. > Эн о,сн,он 


Э) 2 
eon = 25.48310х0н, он) + 4.2821(х0н, он) 1.6354х0н.он (2.7) 


EAC = 104.9956(229 on) + 20.9052(x85 on) 


+2.6349х0н он + 0.4601 (Т2.8) 


Total mass balance 


Total enthalpy flux density 


Conductive heat flow density 


Electrical potential gradient 


Activities (Flory—Huggins 
activity model) 


Conversion of concentrations 


Phase equilibrium (CC-CCP) 


Pressures within membrane are not discussed! 
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Table 3 


Model equations and parameters of models І, Па/ and Ш (both catalyst layers) 


Components 


Model I: anode catalyst layer (AC) 


Model I/model II: cathode catalyst layer (CC) 


Index j =H 20, CH30H, СО» (pure liquid phase) 


Index j= №, Оз, H20, СО» (pure gas phase) 


5 ч Perry F = F ртс 

Reaction mechanism CH30H + H20 < CO? + 6H™ + бе (T3.1) 1.502 + 6H" + 6e —+3H20 (73.12) 

CH3OH + 0.505 9S CO, + 4Ht + 4e7 (T3.13) 
ас^С nO? | AD_gAD = nl мо ассо nM" м.м =n | cdo 
Component mass balance ар = ТАС АС + Vajra (Т3.2) a „СС ICC + рс је + Veross,j”cross (T3.14) 
pores pores 

Stoichiometric constants Va,CH30H = —1; va, H20 = -1; Va,CO = +1 005 = —1.5; Ve,Ny 0; уно = +3; 
YVe,CO7 = 0; Veross, O) = —0.5; Veross,Ny = 0; 
Veross, H20 = 0; Veross,CO7 = +1 

Charge balance 0 = еп — MI Mig With icon = ia = d^ РАС 6 (13.3) 0 = ice – 1м. ам Withicel = (icross + ic); 


Rate expressions 


E AC a,6F AC (CL = a)6F 
Ta = ka е exp (= Па — Xe, ©XP | — prac "a 


withag = 0.1, ka =6 x 1073 тої тт 57! (73.4) 


ic = 0266,06; icross = de Ses 4 Freross (73.15) 
cc 1.5 
k Po, осбЕ (1 —ac)6F 
To = ke exp | — Ne — ех oe 
л 10° pa RTCC ° £ RTOC "° 


ke=2.7 x 1072 mol m~3 s 


M 
"CH, 0H|M=aM 


cross = (Т3.17) 
dee os 
Components Model I/model Ш: anode catalyst layer (AC) Model III: cathode catalyst layer (CC) 
Index j = Н20, CH30H, СО» (pure liquid phase) Index j=N2, O2, H20, СО» (pure gas phase) 
T, To 
Reaction mechanism CH30H + prs pico +4Ht + 4е7 (Т3.5) 0.505 + рро (Т3.18) 
Ru,r, = fre 
PtCO + НО 2 2С05 +2H* + 2e7 + Pt (Т3.6) PtO +2Н? + 2е7 "ре + НОУ (Т3.19) 
ге 
CH30H + P- PICO +4H* + 4е7 (T3.20) 
To, 
рю + РСО- 2р + CO2 (13.21) 
4 
dcAC пАР| ap рар = | Moal, jřal + Yar, jra2 асСС „Ммм =n | ср 
$ j z а [тед j j z" =d j 2-2=0 
Component mass balance = (T3.7) Vek, jřck (13.22) 
dt АС ДАС dt РАС, gcc k 
pores pores 
k=1 
doS = doce x 
СО _ falta ac _ =2 СО _ 103—704 
КТ, WAC with Not 0.11 molm (13.8) а = KE (T3.23) 
pt pt 
dog Гор +Te2 — Tea `СС 2 
= = with № = 0.11 тот (T3.24) 
dt МСС pt 
р 
Stoichiometric constants Val,CH30H = —1; val, H20 = 0; val,CO, = 0; all vek,j = O except for : 1,0) = —0.5; Ve2,H)0 = =1; Ve4,CO7 = +1 
Va2,CH30H = 0; Ya2,Hy0 =—1; %02,С0) = +1 
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account for the influence of the surface coverage of the platinum 
catalyst with adsorbed CO (Langmuir adsorption/desorption 
kinetics). 

The cathode in model II is identical to that in model I. 


3.3. Model III: extended cathode kinetics with adsorbed 
CO and O 


On the cathode platinum catalyst, very similar processes as 
on the anode side take place with respect to the methanol, which 
crosses over through the PEM. Therefore it seems appropriate to 
use similar reaction schemes and rate expressions like the ones 
formulated in model II for the anode side. 

It is assumed, that oxygen and methanol compete for the free 
platinum catalyst surface (Table 3, lower right quarter): oxygen 
is adsorbed in a non-electrochemical reaction forming Pt-O (Eq. 
13.18), methanol is adsorbed forming Pt-CO, as well as protons 
and electrons (Eq. T3.19). Both steps are assumed to be fast. The 
adsorbed oxygen (Pt-O) can then take up protons and electrons 
to form free platinum sites and water (Eq. T3.20). In parallel, 
Pt-CO reacts with Pt-O to carbon dioxide and free platinum sites 
(Eq. T3.21). 

Based on this reaction scheme, both surface species (Pt-O 
and Pt-CO) have to be balanced (OS°, Eq. T3.23 and OSS, 
Eq. T3.24), and the contributions of both electrochemical reac- 
tions are considered in the charge balance (Eq. T3.25). The rate 
expression of the first reaction (Eq. T3.26) is of Langmuir type 
(non-electrochemical irreversible adsorption), the electrochemi- 
cal oxygen reduction is formulated as combined Langmuir-Tafel 
kinetics (Eq. T3.27, analogue to the anode side). The methanol 
adsorption is assumed to be an immediate reaction, therefore 
its rate (Eq. T3.28) is proportional to the methanol crossover 
flux density (as in the previous models). Finally, the release of 
carbon dioxide is described by an irreversible Temkin—Frumkin 
desorption kinetics (Eq. T3.29). 

The anode in model III is the same as in model II. 


4. Simulation results 


In this section, in a first step the presented model variants are 
compared to the experimental results. In a second step, model 
Ш is examined more closely for varying operating conditions 
(in terms of cell temperature and cell current density). 


4.1. Comparison of model variants 


For the operating conditions given in Section 2 and a cell 
current density of 56 mA стт 2, the model variants yield cell 
voltage predictions as presented in Fig. 5. It can be seen that 
model I predicts the cell breakdown much earlier than what was 
found in the experiments, while all other models predict the 
sustaining time correctly. Model II can also predict cell voltage 
overshoots, if the cathode charge transfer coefficient is set to the 
relatively high value of a, = 0.905 (model IIb). Model Ш finally 
can predict overshoots with charge transfer coefficients of 0.5 
on both anode and cathode side. All kinetic parameters used in 
the model variants are given in Table 3. 
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Fig. 5. Comparison of simulated and experimental cell voltages (operating con- 
ditions as in Section 2, cell current density ice = 56 тА cm7?). 


The high value for the cathode charge transfer coefficient in 
model IIb is a strong hint that in the cathode reaction, adsorption 
processes play a major role, justifying the formulation of a more 
complex cathode kinetics like that in model III. 

A more detailed insight in the underlying physico-chemical 
processes can be obtained by a closer look at the state vari- 
ables and the methanol crossover flux densities of the model 
variants as presented in Fig. 6 (anode side) and Fig. 7 (cath- 
ode side). In Fig. 6, the anode catalyst layer state variables 
(methanol concentration сең, он» Surface coverage on the Pt cat- 


alyst @AS, overpotential na) and the methanol crossover flux 
density N¢y,04!-M=o at the interface between anode catalyst 


layer and PEM (zM = 0) are plotted over time. 
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Fig. 6. Comparison of simulation results (anode side) of model variants (oper- 
ating conditions as in Section 2, cell current density ісе = 56 mA cm7?). 
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Fig. 7. Comparison of simulation results (cathode side) of model variants (oper- 
ating conditions as in Section 2, cell current density ice =56mA cm7?). 
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One can see that the variables related to (methanol) mass 
transport (сан,он апа пон) yield very similar results for all 
three models (Fig. 6). This could be expected, as the model vari- 
ants only differ with regard to the assumed electrode kinetics. 
All models predict a fast decrease of the methanol concentra- 
tion in the anode catalyst layer (сё он), slightly (but almost 
negligibly) slowed down by a small reflux of methanol from 
the membrane back into the anode catalyst layer (npon > 0: 


methanol flux towards cathode, пМ он < 0: methanol flux 
towards anode). Less than 20s after the feed switch in all cases 
no more unreacted methanol is available as reactant for the anode 
reactions. 

As model I has no other reactant storages, its anode 
overpotential then sharply increases leading to an immediate 
breakdown of the cell voltage (Fig. 6). In the other model vari- 
ants (ПЛП), the cell can be operated beyond the point where the 
unreacted methanol in the anode catalyst layer is depleted, as 
carbon monoxide (CO) adsorbed to the platinum catalyst forms 
another reactant reservoir. Only when this storage (described 
by the CO surface coverage os) is empty as well, the anode 
overpotential finally increases sharply. Moreover, between the 
feed switch and the final cell breakdown the anode overpoten- 
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Fig. 8. Comparison of simulation results of model Ш to experimental data (operating conditions as in Section 2, except for cell current density ice and anode feed 
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Fig. 9. Comparison of simulation results of model Ш (anode side) (operating conditions as in Section 2, except for cell current density ice and anode feed 
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tial in models II and III also increases much slower than in 
model I. 

Fig. 7 presents the cathode catalyst layer state variables 
and methanol crossover flux density at the interface between 
PEM and cathode catalyst layer (zM =d™) of the model vari- 
ants. Like on the anode side, also here the models do not 
differ significantly in terms of mass transport. A few seconds 
after the feed switch, the methanol crossover flux reaching 
the cathode catalyst (поні 7м_ам) begins to decrease, and 
with it the absolute value of the cathode overpotential. Roughly 
20s after the feed switch (i.e. when the methanol concen- 
tration in the anode catalyst layer, сОҢ,Он» has dropped to 
zero), the methanol crossover flux has ceased completely, and 
the cathode overpotential reaches a new steady state value. 
The decreasing methanol crossover also leads to a parallel 
decrease of the cathode CO coverage (OSS) in model Ш, 
while the oxygen coverage of the cathode catalyst (OSS) is 
increasing accordingly, taking over the sites freed from CO. 
Obviously, all model variants’ cathode processes react with the 
same time constant, they differ only in the steady state over- 
potentials. Therefore, for a simple model, the original cathode 
model with lumped Butler—Volmer type rate expression (Table 3, 


T3.12-T3.17) with appropriately fitted parameters (model Ib) 
is sufficient. 

It is obvious from Figs. 6 and 7, that in order to achieve 
a cell voltage overshoot after the feed switch, the slope of the 
cathode overpotential has to be larger than the slope of the anode 
overpotential. As soon as the cathode arrives at its new steady 
state, the cell voltage reaches its maximum, and then declines 
due to the continuing increase of the anode overpotential. 

It can be seen that models ПЬ and Ш predict the maximum 
of the overshoot earlier than what is observed in the experi- 
ments (Fig. 5). The result would be better if the dead time fgeaa 
accounted for in the experiments would be larger than the esti- 
mated 20s. In fact, this estimation might be too low, as in the 
first seconds after the feed switch the experimental setup is in 
a somewhat transitional state, slowly readjusting anode pres- 
sure and flow rate, while їдеаа was estimated assuming constant 
flow rate without interruptions (see Section 2). The measured 
cell voltage curves (Figs. 1—3) could also support a somewhat 
longer dead time fyeaq Of up to 25—305. 

Nonetheless, such an additional experimental time delay can- 
not fully explain the deviations between times at which the 
maxima of the overshoots are reached in the simulations and the 
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Fig. 10. Comparison of simulation results of model III (cathode side) (operating conditions as in Section 2, except for cell current density ice, and anode feed 


temperature Т^). 


experiments. On one hand, the time constant of the methanol 
crossover might be too low in the models: in this case the trans- 
port coefficients for methanol in the PEM need to be adjusted. On 
the other hand, the cathode overpotential might be less dependent 
on the methanol crossover flux: this could be accounted for by 
adjusting the cathode rate constants, and/or the rate expression 
for the electrochemical step (7,2, Table 3, T3.27). 


4.2. Analysis of model III for different operating conditions 


Figs. 8-10 present simulation results from model Ш for 
two different anode feed temperatures (Т^Ё =60, 75°C) and 
two different current densities (ice = 56, and 124 mA cm7?). 
Fig. 8 compares cell voltages and overpotentials, Fig. 9 com- 
pares the most important anode state variables and mass fluxes, 
while Fig. 10 does the same for the cathode side. The following 
observations and resulting conclusions can be summarized. 


4.2.1. Observation (1) 

The time till the cell breaks down is predicted almost cor- 
rectly by the model (Fig. 8). This is a proof, that the sum of 
internal reactant storages (which are responsible for the sus- 
tained operability after the cut off supply of fresh methanol) is 


quite realistic. Nonetheless, this observation does not give any 
hints on the type and quality of different reactant storages. 


4.2.2. Observation (2) 

The breakdown of the cell voltage is predicted more steeply 
than what was found in the experiments (Fig. 8). In the simu- 
lations the final breakdown slopes of the cell voltage are given 
by the slopes of the respective anode overpotentials. Therefore, 
the model descriptions of the underlying phenomena (transport 
of methanol to the anode as well as anode reaction mechanism 
and reaction rates, Fig. 9) might need some further refinement. 


4.2.3. Observation (3) 

The maxima of the overshoots are predicted too early, and 
the relative overshoots w.r.t. the initial steady state are predicted 
higher than what was found in the experiments (Fig. 8). In the 
simulation, the different behaviour of both the anode and cath- 
ode overpotentials (i.e. their time constants) is responsible for 
the overshoots, as can be seen from Fig. 8. The maximum cell 
voltage is reached when the cathode overpotential has reached 
its new steady state due to the ceasing methanol crossover. Obvi- 
ously, if the methanol crossover (nty,0H|ZM=aM in Fig. 10) 
faded more slowly, the maximum of the overshoot would appear 
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later, and presumably be also less pronounced. Therefore it 
might be beneficial to check the transport properties of the model 
(in the PEM and in the anode diffusion layer), as well as whether 
the experimental dead time is correct (see above sections). If in 
the model the transport of methanol within the anode diffusion 
layer (AD) should be slower due to accounting for the block- 
age of pores by carbon dioxide gas, cCH,OH would decrease 
slower leading to a slower decrease in the methanol crossover 
and thus also in the absolute value of the cathode overpotential. 
The same would be true if the methanol transport inside the PEM 
was slightly slower. As the neglection of a carbon dioxide gas 
phase in the anode pore structures is a significant model simpli- 
fication compared to the real situation, it can be assumed that 
this is the main reason for the presumably too fast decrease of 
the methanol concentration in the anode catalyst layer. In any 
case, if model adjustments on the mass transport side should be 
necessary, the internal reactant storages (especially the number 
of adsorption sites on the anode catalyst, Nec) would have to be 
slightly reduced to achieve the same overall breakdown times as 
in the present simulations. 


4.2.4. Observation (4) 

The model predicts cell voltage overshoots for all four cases 
shown in Figs. 8—10, but in the experiments for the higher 
cell current density (ice = 124 mA cm7?) no overshoots were 
observed. Nonetheless, the model predicts the relative height 
of the overshoots to decrease with increasing current density. 
Adjusting the rate constants can decrease the overshoots in gen- 
eral, therefore a disappearance for high current densities seems 
possible. To correctly predict the experimental behaviour for 
the higher current density, the absolutes of the slopes of the 
cathode and anode overpotentials, respectively, after the feed 
switch would have to be identical. Then the cell voltage would 
remain constant till the cathode overpotential approaches its new 
steady state, while the anode overpotential would still increase, 
thus leading to a breakdown of the cell voltage in the observed 
manner. 


4.2.5. Observation (5) 

Finally, from the steady state simulation results prior to the 
feed switch, one can see that the model predicts an increase 
of the cell voltage with higher cell temperature, while in the 
experiments a contrary behaviour is observed (Fig. 8). In the 
presented model, the rate constants are fitted for an anode 
feed temperature of 60°C, no additional (Arrhenius) term 
accounts for changes in the reaction rate constants with tem- 
perature (see Table 3, rate Eqs. T3.4, T3.10, T3.11, T3.16 and 
T3.26-T3.29). A modification of the rate constants by intro- 
ducing such an Arrhenius term will be part of the next steps in 
model refinement. For the anode side, the additional parameters 
(activation energies) have been determined from experiments 
recently [8]. 


5. Conclusions 


The dynamic response of the DMFC to step changes in 
the anode feed methanol concentration was analysed experi- 


mentally, and a dynamic process model was formulated. The 
experiments revealed, that for switching the feed from methanol 
solutions to pure water the cell voltage breaks down later than 
expected, and for low current densities even shows an overshoot 
behaviour. 

The here presented model shows good quantitative agree- 
ment to steady state experiments, and also a good qualitative 
agreement to the aforementioned dynamic step-down exper- 
iments (long sustaining times and overshoots of the cell 
voltage). 

According to the model analysis, the overshoot phenomenon 
can be explained by different response times (time constants) 
of the anode and cathode overpotentials, respectively. After 
rinsing the anode compartment, the methanol concentration 
in the anode catalyst layer drops rapidly due to the ongoing 
methanol consumption by the electrochemical reactions and 
methanol transport towards the anode compartment. Conse- 
quently, the methanol crossover to the cathode decreases, leading 
to a decrease of the absolute cathode overpotential. The anode 
overpotential increases quite slowly, and the cell voltage does 
not break down when the reactant methanol is completely con- 
sumed. The latter is due to reaction intermediates of the anodic 
methanol oxidation adsorbed to the anode catalyst (e.g. CO), 
which form a second reactant reservoir. This enables a prolonged 
operation of the cell for a few seconds up toa few 10 s (depending 
on current density). 

The reason for the observed overshoots therefore seems to 
be, that the absolute cathode overpotential (mainly determined 
by methanol crossover) decreases faster, than the anode overpo- 
tential increases in this transient situation. 

The analysis has also revealed several hints for further model 
improvement: 


e More realistic anode diffusion layer mass transport account- 
ing for the influence of carbon dioxide gas bubbles blocking 
part of the available pores. 

e Validation of the experimental dead time. 

e Adjustment of reaction rate expressions and corresponding 
parameters. 


On the other hand it can be concluded, that for a reduced 
model it might be sufficient to use a simplified cathode reaction 
mechanism (as in model I/II) with Butler—Volmer or even Tafel 
type rate expressions. Nonetheless, a realistic anode reaction 
mechanism (e.g. models П/Ш), seems to be necessary in order 
to correctly predict the dynamic DMFC behaviour. 
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